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CABINET PLASSERAUD 


PROCEDE ET DISPOSITIF ^IDENTIFICATION ADAPT AT1VE. ET 
ANNULEUR P'ECHO ADAPTATIF 
METTANT EN CEUVRE UN TEL PROCEDE 


La presente invention concerne de fa$on generate les techniques 
d' identification adaptative. Elle trouve une application particuliere, non 
exclusive, dans les annuleurs d'echo utilises en telecommunications. 

identification adaptative de systemes lineaires caracterises par leur 
reponse impulsionnelle a ete largement etudiee, et de nombreuses solutions 
algorithmiques ont ete proposees dans la litterature. 

On considere le probleme general de Identification directe par filtrage 
adaptatif transverse, ce qui correspond a la grande majorite des applications 
pratiques. 

La figure 1 montre un systeme a identifier 10 auquel est applique un 
signal x< variant dans le temps. On note z< la reponse du systeme 10 au signal 
d'entree La mesure de la reponse s'accompagne inevitablement de 
I'adjonction dune composante de perturbation b t appelee bruit d'observation. Ce 
bruit d'observation b t peut comprendre du bruit a proprement parler (bruit blanc 
ou bruit routier par exemple), mais egalement du signal utile. La composante b t 
est appelee bruit d'observation dans la mesure ou elle perturbe I'otpervation de 
la reponse Zj. L'additionneur 12 symbolise sur la figure Tadjonction de la 
composante de perturbation b t , supposee additive, a la reponse Le signal 
d'observation mesure y t est ainsi la reponse tfun systeme reel 14 incluant le 

systeme a identifier 10 et I'additionneur 12. 

Le dispositif ^identification adaptative 16 re?oit le signal d'entree ^ sur 

une premiere entree E1 , et le signal d'observation y t sur une seconde entree E2. 
Les signaux x t et y t sont amplifies, filtres et numerises a I'entree du dispositif 16 
par des elements classiques non representes. Le dispositif ^identification 
adaptative 16 comporte un filtre d'identification 18 consistant en un filtre 
programmable a reponse impulsionnelle finie ou FIR (« Finite Impulse 
Response ») notee Hj^ = hj_ r - hj^ 1 ), ou (.) T designe la transposition 
matricielle. Les coefficients du filtre d' identification 18 sont adaptes pour que 
cette reponse impulsionnelle hJ\, soit representative de la reponse 


impulsionnelle du systeme a identifier 10. Le filtre 18 regoit le signal tfentree 
numerise Xj, et delivre une estimation z t de la reponse du systeme 10. 

Un soustracteur 20 retranche cette estimation z t du signal 

^observation numerise y t pour fournir un signal d'erreur e t . Ce signal d'erreur e t 

peut etre vu comme une estimation de la composante de perturbation b t . 

Une unite 22 de mise a jour du filtre ^identification adapte les 
coefficients du filtre 1 8 sur la base du signal d'entree ^ et du signal d'erreur e t , en 

prenant en compte generalement un pas ^adaptation ji. 

De nombreux algorithmes ont ete proposes pour determiner de fagon 
automatique les coefficients du filtre adaptatif 18. La mise en oeuvre de ces 
dispositifs amene generalement leurs concepteurs a etudier le compromis 
devant etre realise entre rapidite de convergence de Talgorithme, simplicite de 
controle et de mise en ceuvre, complexite arithmetique, et stabilite numerique. 

L'algorithme LMS (« Least Mean Square ») est Talgorithme le plus 
repandu pour adapter de fagon continue au cours du temps la reponse 
impulsionnelle d'un filtre ^identification de type FIR. Un tel algorithme permet 
d'implementer de maniere particulierement efficace le filtre de Wiener a L 
coefficients qui minimise, selon Tapproximation stochastique, la valeur 
moyenne de la puissance de I'erreur de filtrage. II est defini par les relations : 

*t = y t -x t T H t _ 1 (1) 

H t =H t .,+ l ie t X t (2) 
ou ^ = (xt,x t -i,--- 1 xt»L+i) T represente le vecteur des L derniers echantillons 
du signal d'entree, et n represente le pas d'adaptation de ralgorithme. Les 
principaux atouts de cet algorithme sont sa faible complexite numerique, sa 
simplicite d'implementation et sa robustesse aux erreurs. Malheureusement, 
lorsque des signaux fortement correles (tels que des signaux de parole) sont 
utilises pour exciter le systeme inconnu, cet algorithme possede une vitesse de 
convergence qui se degrade rapidement. 

Pour s'affranchir de ces inconvenients, il est souvent fait usage 
d'une version specifique de Talgorithme LMS integrant un pas d'adaptation 
parametrise. Cet algorithme correspond alors a une version normalisee du 
LMS, ou NLMS (« Normalized Least Mean Square ») dans laquelle les 
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coefficients du filtre adaptatif sont mis a jour selon la relation suivante : 


un filtre de type FIR d'ordre inferieur ou egal a L la relation (3) peut s'ecrire : 


ou AH t =H opt -H t represente Terreur dans les coefficients estimes du filtre a 

Iteration t. Cette expression correspond a une interpretation geometrique de 
ralgorithme NLMS. Dans les cas ou n*1 , r equation (4) correspond a une 

projection relachee du vecteur AH t _i sur un sous-espace affine completement 
determine par la matrice comprise entre les crochets et par la connaissance du 
decalage a I'origine donne par le dernier terme de la relation (4). 

Afin d'apporter de nouveaux algorithmes offrant une rapidite de 
convergence superieure a celle de ralgorithme NLMS, plusieurs approches ont 
ete proposees dans la litterature (implementation dans le domaine frequentiel, 
filtrage en sous-bandes...)- On considere ci-dessous celles basees sur une 
modification de la direction de projection du NLMS, ainsi que celles basees sur 
I'utilisation d'un filtre blanchissant. 

Amelioration de la convergence oar modification de la direction de projection 

II est possible d'ameliorer la convergence du NLMS en modifiant la 
direction de la projection mentionnee ci-dessus. Cette analyse est a Torigine 
des algorithmes de projection affine, ou APA (« Affine Projection Algorithm ») 
qui sont bases sur une projection d'ordre multiple egal a P. Ces algorithmes 
possedent de ce fait de bien meilleures proprietes de convergence sur des 
signaux correles en comparaison avec ralgorithme NLMS (qui correspond au 
cas limite ou P=1). Les algorithmes APA d'ordre P (voir K. Oseki et al., « An 
adaptative algorithm using an orthogonal projection to an affine subspace and 
its properties » t Electronics and communications in Japan, 1984, Vol. 67-A, 
n° 5, pages 19-27) sont caracterises par la mise a jour des coefficients du filtre 
d'identification 18 suivant les relations : 



En supposant que le filtre optimal H opt du systeme inconnu 10 soit 



(4) 


e t f p = Y t,p - 2St,p H t-1 


(5) 



Jpfetp^jp)" 


(6) 


(7) 


(8) 


represente la matrice inverse generalisee de Moore-Penrose d'ordre LxP. Ces 
equations de mise a jour des coefficients du filtre ^identification permettent de 
montrer que le vecteur d'erreurs d 1 estimations a posteriori e^^f* est egal a : 


En considerant le pas d'adaptation de ralgorithme egal a I'unite, 
1'algorithme de projection affine d'ordre P annule les P erreurs a posteriori 
definies a la relation (9). Cette derniere propriete explique le tres bon 
comportement en vitesse de convergence de cet algorithnne. 
Malheureusement, dans sa version de base decrite par les relations (5) et (6), 

la complexity theorique de tels algorithmes est de I'ordre de 2LP + K jnv P 2 ou 
Kj n v represente une constante associee au calcul de la matrice inverse 
intervenant dans I'expression (6), les parametres L et P designant 
respectivement le nombre de coefficients du filtre ^identification et I'ordre de 
projection. 

Afin de reduire cette complexite initiale, plusieurs versions rapides de 
ces algorithmes ont ete proposees en segmentant la pseudo-matrice 
d'autocorrelation de fagon similaire aux algorithmes des moindres carres 
recursifs rapides. De telles techniques permettent ainsi de reduire la 
complexite initiale a des valeurs plus raisonnables de Tordre de 2L + 20P (voir : 
Steven L. Gay, « A fast converging, low complexity adaptative filtering 
algorithm », Proceeding of the 3 rd International Workshop on Acoustic Echo 
and Noise Control, Plestin-Les-Greves, France 1993, pages 223-226 ; Steven 
L. Gay, « Fast projection algorithms with application to voice echo 
cancellation », Ph.D. Dissertation of the State University of New Jersey, USA, 
1994 ; Steven L. Gay et al., « The fast affine projection algorithm », Proceeding 
of ICASSP'95, pages 3023-3026, 1995; M. Montazeri, « Une famille 
d'algorithmes adaptatifs comprenant les algorithmes NLMS et RLS : application 
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a I'annulation d'echo acoustique », These de Doctorat de I'Universite de Paris 
Sud t 1994 ; M. Tanaka et al., « Reduction of computation for high-order 
projection algorithm », Electronics Information Communication Society Autumn 
Seminar, Tokyo, 1993). 

5 Amelioration de la convergence oar insertion d'un filt raoe blanchissant 

Plusieurs travaux de recherche ont ete consacres a Tetude de 
Amelioration des performances des systemes d' identification adaptative a 
travers I'emploi de structures predictives (voir : M. Mboup et al., « LMS coupled 
adaptive prediction and system identification : a statistical model and transient 

10 analysis », IEEE Transactions on signal processing, Vol. 42, n° 10, octobre 
1994, pages 2607-2615 ; S. Benjebara, « Caracteristiques des signaux et 
capacite de poursuite des non-stationnarites aleatoires : apport des schemas 
predictifs et multiresolutions », These de I'Universite des Sciences, des 
Techniques et de Medecine de Tunis II, Tunis, 1997). Ces etudes ont permis 

15 de degager principalement deux structures (symetrique ou non-symetrique) de 
pre-blanchiment du signal d'excitation du filtre par une technique de filtrage 
adaptatif. Le principe general de traitement de la structure non-symetrique est 
schematise a la figure 2. 

y ne telle structure repose essentiellement sur une demarche 

20 empirique visant a modifier le signal utilise dans la mise a jour des coefficients 
du filtre adaptatif et a le transformer de fagon a reduire le conditionnement de 
sa matrice ^autocorrelation (rapport entre les valeurs propres maxi males et 
minimales de la matrice d'autocorrelation de ce signal). De ce fait, la mise a 
jour des coefficients est realisee par le module d'adaptation 22 a partir du 

25 signal disponible en sortie d'un circuit de prediction lineaire 24, la prediction 
etant realisee a partir de M coefficients. Ualgorithme applique par le module 22 
pour mettre a jour les L coefficients du filtre d' identification 18 correspond au 
LMS (relation (2)) ou au NLMS (relation (3)). De meme, le circuit de prediction 
24 est implements sous la forme d'un filtre adaptatif a M coefficients dont la 

30 mise a jour est realisee a Taide d'un algorithme LMS. 

Les etudes consacrees a Tanalyse des performances obtenues a partir 
de ce type de structure ont notamment mis en evidence Texistence tfune 
interaction tres forte entre les modules de prediction et d'adaptation du 
systeme. Ceci entralne notamment une interdependance entre les equations 
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cTadaptation des coefficients du predicteur adaptatif 24 et celles du filtre 
d'identification 18. Cette forte dependance se retrouve egalement dans le choix 
des deux pas ^adaptation |i p et p. H) ce qui engendre une zone de stabilite tres 

reduite de la structure complete. 
5 En consequence, les auteurs mentionnent une instability du systeme 

complet d 1 identification par structure predictive, ce qui limite les ordres de 
prediction utilises dans les realisations a des valeurs toujours inferieures a 
quatre afin d'assurer une relative stabilite, ceci limite les performances de 
identification : on est amene a choisir des pas d'adaptation up e * Uh tlf ss 

10 faibles, non compatibles avec Tobjectif ^amelioration de la rapidite de 
convergence. 

L'examen des solutions proposees dans la litterature pour ameliorer la 
rapidite de convergence des algorithmes d' identification adaptative conduit a 
conclure d'une part que les algorithmes bases sur une modification de la 

15 direction de projection restent complexes en termes de nombre d'operations 
arithmetiques requises pour la mise a jour des coefficients du filtre lorsque 
I'ordre de projection P est eleve, et d'autre part que les schemas a base 
^identification par structure predictive restent tres delicats a controler et offrent 
un gain de convergence reduit du fait des faibles ordres de prediction M utilises 

20 en pratique afin d'assurer une relative stabilite de la structure globale. 

La presente invention a pour but de proposer un procede 
d'identification adaptative qui presente de bonnes proprietes de convergence et 
qui soit d'une mise en ceuvre relativement simple et d'une complexite 
arithmetique limitee. 

25 Uinvention propose ainsi un procede d'identification adaptative pour 

estimer la reponse d'un systeme a un signal tfentree, dans lequel : on re?oit 
d'une part le signal d'entree et d'autre part un signal d'observation dont une 
composante est ladite reponse au signal d'entree ; on determine un signal 

d'erreur e t a un instant t selon Texpression e t = y t - X J H t-1 , ou y t designe la 
30 valeur du signal d'observation a I'instant t, H t _^ est un vecteur-colonne 
compose de L coefficients d'un filtre d'identification a reponse impulsionnelle 
finie representative de la reponse impulsionnelle du systeme, et 

X^= (xt,x t _ 1 ,---,x t _ L+1 ) est un vecteur-ligne compose des valeurs x t , x t _ 1 ,..., 


x t . L+1 du signal d'entree a I'instant t et aux L-1 instants precedents ; et on 
adapte les L coefficients du filtre d' identification en tenant compte du signal 
d'entree et du signal d'erreur. Suivant I'invention, on obtient des parametres de 
prediction du signal d'entree, minimisant I'energie d'un residu de prediction sur 
des trames successives du signal d'entree, et 1'adaptation des L coefficients du 
filtre ^identification consiste a ajouter au vecteur-colonne H t .-, un vecteur- 

colonne proportionnel a _ et U t , ou U t est un vecteur-colonne compose 

des L valeurs du residu de prediction du signal d'entree a I'instant t et aux L-1 
instants precedents, et X un coefficient positif ou nul. 

Dans des modes d'execution preferes du procede : 

- les trames du signal d'entree ont une duree d'au moins 5 ms ; 

- les trames du signal d'entree presentent des recouvrements mutuels ; 

- on effectue une analyse par prediction lineaire d'ordre P-1 du signal 
d'entree, d'ordre P-1 de preference egal a 5 ou plus, et le vecteur- 

P-1 

colonne U t a pour expression U t =X t + ^a q X t _ q , ou les a q designent 

q=1 

les coefficients de prediction lineaire resultant de ladite analyse. L'ordre 
P_1 de | a prediction lineaire peut egalement etre une fonction d'une 
.stationnarite estimee du signal d'entree ; 

- le signal d'entree est un signal de parole restitue par un decodeur a partir 
d'un flux binaire d'entree, et les parametres de prediction du signal 
d'entree sont extraits du flux binaire d'entree par le decodeur qui peut en 
outre fournir le residu de prediction pour Tadaptation des L coefficients du 
filtre ^identification. 

Un autre aspect de I'invention se rapporte a un dispositif ^identification 
adaptative d'un systeme auquel est applique un signal d'entree, comprenant : 
une premiere entree recevant le signal d'entree ; une seconde entree recevant 
un signal d'observation dont une composante est une reponse du systeme au 
signal d'entree; un filtre ^identification a reponse impulsionnelle finie 
representative de la reponse impulsionnelle du systeme ; un soustracteur 
produisant un signal d'erreur e t ayant I'expression (1) ci-dessus ; des moyens 
pour obtenir des parametres de prediction du signal d'entree, minimisant 
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I'energie d'un residu de prediction sur des trames successives du signal 
d'entree ; et des moyens d'adaptation mettant a jour les L coefficients du filtre 
d'identification en ajoutant au vecteur-colonne H t . 1 le vecteur-colonne 


proportionnel a 


x t T u t +x *' 

Un autre aspect de I'invention se rapporte a un annuleur d'echo 
adaptatif pour enlever d'un signal de retour des composantes d'echo d'un 
signal direct, comprenant un dispositif d'identification du type enonce ci- 
dessus, dont la premiere entree regoit le signal direct en tant que signal 
d'entree et dont la seconde entree resoit le signal de retour en tant que signal 
d'observation, le signal d'enreur constituant le signal de sortie de I'annuleur 
d'echo. 

D'autres parti cularites et avantages de la presente invention 
apparaitront dans la description ci-apres d'exemples de realisation non 
limitatifs, en reference aux dessins annexes, dans lesquels : 

- les figures 1 et 2, precedemment commentees, sont des schemas 
synoptiques de dispositifs d'identification adaptative de Tart anterieur ; 

- la figure 3 est un schema synoptique d'un annuleur d'echo incorporant un 
dispositif d'identification adaptative selon I'invention ; et 

- la figure 4 est un schema synoptique d'une autre forme de realisation 
d'un annuleur d'echo selon I'invention. 

Dans le cas de signaux stationnaires, on peut etablir la relation 
suivante entre les coefficients a 1t ... f a p ^ de prediction lineaire d'ordre P-1 du 

signal d'entree et I'inverse R~p de la matrice d'autocorrelation 
R t p = *\pXjp du signal (la matrice Xj p etant definie par la relation (7)) : 


(10) 


( 1 ^ 


{^\ 

0 


= R t,p 





avec 


p-1 

U t = X t + £>q X t-q 
q=1 


(11) 


Si on considere I'algorithme APA d'ordre P defini plus haut avec un pas 
d'adaptation n=1, la nullite des erreurs a posteriori (equation (9)) simplifie le 


vecteur d'erreurs a priori (5) : e t p = (e t( 0,- I 0) T l de sorte que la formule (6) 
d'adaptation du filtre d'identification devient, compte tenu de la relation (10) : 


H t =H M + U t (12) 

X^Ut 

Sur cette base, on peut definir une variante de I'algorithme APA d'ordre 
P, ci-apres appelee « algorithme pseudo-APA d'ordre P », dans laquelle on 
reintroduit un pas d'adaptation \x et eventuellement une constante de 
regularisation positive X. Dans I'algorithme pseudo-APA, le signal d'erreur est 
defini par la relation (1), et le filtre d'identification est mis a jour selon : 

H t =H M + U t (13) 

X t T U t +X 

Le pas d'adaptation \x est typiquement compris entre 0,1 et 1, et X > 0. 
Par rapport a la structure d'identification par prediction decrite precedemment, 
I'algorithme pseudo-APA d'ordre P presente deux differences irnportantes : 

- les coefficients de prediction a q sont evalues non plus au rythme des 

echantillons mais a un rythme de blocs, lis correspondent a la 
minimisation de Tenergie de prediction lineaire sur des trames de taille N 
echantillons. Des techniques efficaces d'implementation peuvent etre 
utilisees afin d'estimer I'ensemble {a q ; q=1, P-1} des coefficients du 

filtre predicteur. 

- Tequation (13) de mise a jour des coefficients du filtre d'identification 
Ht.! ne correspond pas a celle utilisee dans les travaux anterieurs. En 
effet, ces travaux reposent essentiellement sur une approche empirique 
visant a pre-blanchir le signal d'excitation afin d'ameliorer la vitesse de 
convergence des modes, celle-ci etant proportionnelle au 
conditionnement de la matrice d'autocorrelation du signal d'excitation. De 
ce fait, les algorithmes LMS ou NLMS proposes par les auteurs pour la 
mise a jour des coefficients du filtre d'identification (et egalement du 
predicteur) ne correspondent pas a celui de la relation (13). Ces 
algorithmes different par le terme de normalisation present dans 
('equation de mise a jour des coefficients du filtre, le terme de la forme 
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[U t U t j du NLMS etant remplace dans la relation (13) par un terme 


Le procede selon I'invention resout la problematique frequemment 
posee par la forte correlation des signaux a traiter. En effet, dans ce dernier 
cas, les algorithmes classiques tendent a perdre leurs bonnes proprietes de 
convergence. Le procede selon I'invention permet intrinsequernent de 
conserver ces bonnes proprietes de convergence meme en presence de 
signaux fortement correles, ce qui est souvent le cas en pratique. II permet, en 
outre, de conserver une complexite arithmetique limitee ce qui represente un 
tres grand avantage pour Timplementation en temps reel de ces algorithmes 
d' identification sur processeurs de traitement du signal. 

L'algorithme pseudo-APA est a la base d'une nouvelle famille de 
dispositifs d'identification adaptative qui peuvent etre utilises dans divers 
domaines d'application (annulation d'echo, egalisation de canaux de 
propagation, commande automatique de processus, etc.). On considere ci- 
apres, sans que cela soit limitatif, I'exemple de I'annulation d'echo acoustique. 

La figure 3 montre un annuleur d'echo adaptatif incorporant un dispositif 
d 1 identification adaptative 16 selon I'invention. L'annuleur d'echo est associe a 
une installation telephonique mains-libres. Le signal cfentree x l , regu sur I'entree 

E1 du dispositif, est le signal direct destine au haut-parleur 11 de Installation 
mains-libres. Le signal d'observation y t , regu sur I'entree E2 du dispositif, est le 

signal de retour recueilli par le microphone 13 de Installation mains-libres. Ce 
signal d'observation y t inclut des composantes d'echo ^ du signal direct, et des 

composantes de perturbation b t pouvant inclure du bruit et de la parole emise par 

le locuteur. Dans ce cas, le systeme a identifier consiste en le ou les trajets 
d'echo entre le haut-parleur 1 1 et le microphone 13. 

Le signal de sortie de l'annuleur d'echo est le signal d'erreur e t delivre 

par le soustracteur 20 a partir du signal d'observation y t et de Testimation d'echo 

z t produite par le filtre d' identification 1 8. 

Les fenetres temporelles du signal d'entree sont gerees par un 

module 30 qui memorise au moins les L derniers echantillons Xj, Xj.^..., x t _ L+1 . 

Ces L echantillons, qui forment le vecteur X t , sont fournis : (i) au filtre 



d' identification 18 qui precede au filtrage conformement au dernier terme de la 
relation (1 ) ; (ii) au module d'adaptation 22 qui met a jour les coefficients du 
filtre 18 suivant la relation (13) ; et (iii) au filtre de prediction 32 qui fournit au 
module d'adaptation 22 le vecteur U t defini par la relation (11), compose des L 

dernieres valeurs du residu de prediction. 

A chaque instant t, ii suffit que le filtre 32 produise la premiere 
composante du vecteur U t , c'est-a-dire la valeur courante du residu egale a 


apres reception des precedents echantillons. 

Dans I'exemple represents sur la figure 3, le filtre de prediction 32 a 
une structure classique en treillis, les coefficients de prediction a q (1<q<P) 

etant representes par les P-1 coefficients de reflexion assbcies r^..., rp.-,. Les 

coefficients de reflexion r t sont par exemple obtenus par un module 34 

conformement a Talgorithme bien connu de Levinson-Durbin, base sur les 
autocorrelations <J>(i) du signal d'entree calculees par un module" 36. 
L'algorithme de Levinson-Durbin execute par le module 34 a Texpression 
rappelee ci-dessous : 


P-i 

x t + ^ a q x t-q i les autres composantes ayant ete calculees et memorisees 
Q=1 


e(0) = m 

Pour i allant de 1 a P-1 , faire : 



a = -n 


E(i) = (l- ri 2 )E(i-1) 


Pour j allant de 1 a i-1, faire : 



Fin Pour 


Fin Pour 


Les coefficients de reflexion r s fournis au filtre 32 sont ceux obtenus a 
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I'issue de I'iteration P-1, de merne que les coefficients a p de la relation (11) 
(a q = aq" 1 ). La quantite E(P-1) est I'energie de I'erreur residuelle de prediction. 
Les coefficients de correlation <|>(i) sont calcul es par le module 36 

selon : 

N-1 

<K')=2x(j).x(H) (14) 
j=0 

ou xQ) = x t _j.f(j) represente le signal d'entree multiplie par une fonction de 
fenetrage classique telle qu'une fonction rectangulaire, de Hamming ou autre. 

Le calcul des correlations <|>(i) et I'algorithme de Levinson-Durbin, qui 
equivaut a une minimisation de I'energie de prediction E(P-1 ), sont executes 
sur des frames de N echantillons du signal d'entree, N etant un nombre du 
meme ordre de grandeur que la longueur L de la reponse impulsionnelle du 
filtre d'identification 18. 

A titre d'exemple, les signaux etant echantillonnes a une frequence 
Fe = 8 kHz, le filtre 18 peut avoir L = 256 coefficients, et la taille des trames 
peut etre N = 160, c'est-a-dire que les trames sont de 20 ms. Dans une 
application de phonie, on peut considerer que, sur une telle duree, les signaux 
de parole sont quasi-stationnaires, ce qui valide I'une des hypotheses faites 
dans la derivation de I'algorithme pseudo-APA a partir de I'algorithme APA. En 
general, cette duree du trame pourra etre prise superieure a 5 ms. 

Les trames du signal d'entree traitees par les modules d'analyse 
predictive 34, 36 sont de preference recouvrantes, ce qui permet de tenir 
compte du caractere non-stationnaire du signal. L'analyse predictive est alors 
effectuee tous les K echantillons du signal d'entree, avec K < N. A titre 
d'exemple, la duree de recouvrement entre deux trames successives peut etre 
de Pordre de 15 ms. Avec I'exemple numerique precedent, cela correspond a 
K = 40, chaque frame etant compose de quatre blocs successifs. 

Quand N est un multiple de K, ceci permet en outre de simplifier le 
calcul des coefficients de correlation par le module 36. Apres reception de 
chaque bloc de K valeurs du signal d'entree ^ il lui suffit de calculer, pour 
chaque index i, une correlation partielle correspondant aux K termes les plus 
recents de I'equation (14), et de mettre a jour les correlations <|>(i) en leur 
ajoutant les correlations partielles qui viennent d'etre calculees et en leur 
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retranchant les correlations partielles calculees N/K blocs precedemment. 

On notera que d'autres techniques peuvent etre employees pour 
mettre en oeuvre le filtre de prediction. On peut ainsi utiliser une methode autre 
que ralgorithme de Levinson-Durbin pour calculer les coefficients de reflexion rj 

5 ou une autre structure du filtre de prediction basee sur les coefficients de 
reflexion r y On peut egalement utiliser une autre structure connue du filtre de 

prediction basee directement sur les coefficients de prediction a q , ou sur les 

coefficients LAR (« Log-Area-Ratio », LARj = log 10 [(1-rj)/(1+rj)] ), ou encore sur 

les coefficients LSP (« Line Spectrum Pair ») ... 
10 Le dispositif selon r invention offre la possibility de faire varier I'ordre 

P-1 de la prediction en fonction des caracteristiques instantanees du signal 
d'entree. En particulier, on peut evaluer la stationnarite du signal x^ par 

exemple en analysant les coefficients de correlation ou de correlation partielle 
calcules par le module 36, afin d'adopter un ordre de prediction P-1 
15 relativement eleve lorsque le signal d'entree est plutot stationnaire, et un ordre 
de prediction plus faible en presence de non-stationnarites. 



NLMS 

APA exact 

Pseudo-APA 

Complexite theorique 

3L 

2 LP + K,. nv P2 

Max{3L, 

par echantillon (en 
nombre d'operations 



^P + (P'+5).P'} 

elementaires) 




Complexite/echantillon 

768 op. 

> 4096 op. 

824 op. 

avec la configuration : 




Fe = 8kHz 




L = 256 coefficients 

6,14 

> 32,7 Mflops 

6,6 Mflops 

N = 160 (20ms) 

Mflops 



F = P-1 = 8 





Tableau I 


II est possible d'evaluer la complexite arithmetique de ralgorithme 
pseudo-APA. Les caracteristiques pertinentes sont resumees dans le Tableau 
20 I, ou elles sont comparees a celles des algorithmes NLMS et APA exact d'ordre 
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P, dans une configuration typique cle I'annulation d'echo. Dans le tableau, P 
designe I'ordre maximal de Talgorithme pseudo-APA que peut mettre en oeuvre 
le circuit 16, et P' I'ordre effectivement utilise dans I'exemple comparatif. 

On constate que la complexity arithmetique de I'algorithme pseudo- 
5 APA est sensiblement equivalente a celle du NLMS dans le cas d'un filtre 
d' identification comportant L = 256 coefficients, et pour un filtre de prediction 
dont les coefficients sont evalues toutes les 5 ms sur des trames de 20 ms 
(K = 40, N = 160). Linvention propose ainsi des algorithmes approchant les 
performances des algorithmes de projection affine d'ordre eleve (typiquement 

10 superieur a 4) pour une complexity inferieure a celle de I'algorithme exact de 
projection affine d'ordre 2. 

Materiellement, le dispositif ^identification 16 peut etre realise a partir 
de circuits du commerce, notamment des processeurs de traitement de signal 
en temps reel (DSP) a arithmetique flottante couramment utilises (TMS320C3X 

15 et TMS320C4X vendus par la societe Texas Instruments, AD21061 vendu par 
la societe Analog Devices ...). Des processeurs a arithmetique fixe peuvent 
egalement etre utilises, en respectant les precautions usuelles de cadrage 
approprie des donnees dans les traitements. 

Dans la variante de realisation representee schematiquement sur la 

20 figure 4, le signal d'entree Xj est un signal de parole reconstruit par un 

decodeur 40 a partir d'un flux binaire d'entree O. Le flux binaire & a ete produit 
par un codeur dual du decodeur 40, situe dans le reseau de communication ou 
encore dans le terminal utilise par le locuteur distant. 

A titre d'exemple, le codeur/decodeur peut correspondre au codec 

25 RPE-LTP employe dans le systeme GSM de radiotelephonie cellulaire (voir 
Specification GSM 06.01 diffusee par I'European Telecommunications 
Standard institute). Ce codeur, comme une majorite des codeurs audio 
numeriques d'utilisation courante, procede a une analyse par prediction lineaire 
sur des trames du signal de parole correspondant a une zone de quasi- 

30 stationnarite du signal, ce qui repond bien au besoin de la presente invention. 
Dans le cas du GSM, les trames sont de 20 ms. Le codeur inclut dans le flux 
binaire <I> qu'il delivre d'une part des parametres de quantification des 
coefficients LAR evoques plus haut, qui caracterisent le filtre de prediction 
lineaire d'ordre 10 du signal de parole, et d'autre part des parametres de 
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quantification d'un signal d'excitation correspondant au residu de prediction 
lineaire. 

Le decodeur 40 comporte des drcuits qui recuperent le signal 
d'excitation a partir de ses parametres de quantification lus dans le flux <I>, et 
un filtre de synthese inverse du filtre de prediction d'ordre 10, auquel sont 
fournis les coefficients deduits des LAR quantifies. Le filtre de synthese regoit 
le signal d'excitation, et delivre le signal de parole decode x t . Le decodeur 40 

peut done obtenir les parametres de prediction lineaire dans le flux <1>, a la 
place des modules 34 et 36 du dispositif de la figure 3. II peut en outre delivrer 
le residu de prediction (correspondant au signal d'excitation du filtre de 
synthase) formant les composantes du vecteur U t utiles au module 

d'adaptation 22. 

Dans le cadre d'un equipement comportant un decodeur de parole et 
un annuieur d'echo adaptatif, on voit que Tutilisatton de Talgorithme pseudo- 
APA cpnformement au schema de la figure 4 n'introduit pas de complexite 
supplementaire par rapipwl au NLMS, quel que spit I'ordre P. 
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REVENDICATIONS 


Procede d' identification adaptative pour estimer la reponse d'un 


systeme (10) a un signal d'entree, dans lequel : on regoit d'une part le signal 
d'entree et d'autre part un signal d'observation dont une composante est ladite 
5 reponse au signal d'entree ; on determine un signal d'erreur e t a un instant t 

selon I'expression e t = y t -X^H M , ou y t designe la valeur du signal 
d'observation a I'instant t, H t . 1 est un vecteur-colonne compose de L 
coefficients d'un filtre ^identification (18) a reponse impulsionnelle finie 
representative de la reponse impulsionnelle du systeme, et 

10 X^ = (xt,x t _ 1 ,---,x t _ L+1 ) est un vecteur-ligne compose des valeurs x^ x^,..., 
Xt-L+1 du si 9 nal d'entree a Tinstant t et aux L-1 instants precedents ; et on 
adapte les L coefficients du filtre d' identification en tenant compte du signal 
d'entree et du signal d'erreur, caracterise en ce qu'on obtient des parametres 
de prediction du signal d'entree, minimisant I'energie d'un residu de prediction 

15 sur des trames successives du signal d'entree, et en ce que I'adaptation des L 
coefficients du filtre ^identification consiste a ajouter au vecteur-colonne H* * 


compose des L valeurs du residu de prediction du signal d'entree a I'instant t et 
aux L-1 instants precedents, et X un coefficient positif ou nul. 


un vecteur-colonne proportionnel a 



U t , ou U t est un vecteur-colonne 


20 


2. Procede selon la revendication 1, dans lequel lesdites trames du 

signal d'entree ont une duree d'au moins 5 millisecondes. 


3. Procede selon la revendication 1 ou 2, dans lequel lesdits trames du 

signal d'entree presentent des recouvrements mutuels. 


4. Procede selon I'une quelconque des revendications 1 a 3, dans 

25 lequel on effectue une analyse par prediction lineaire d'ordre P-1 du signal 
d'entree, et dans lequel le vecteur-colonne U t a pour expression 


p-1 

U t = X t + X] a qXt-q . ou les a q designent les coefficients de prediction lineaire 
q=l 

resultant de ladite analyse. 

5. Procede selon la revendication 4, dans lequel I'ordre P-1 de la 
prediction lineaire est au moins egal a 5. 

6. Procede selon la revendication 4 ou 5, dans lequel I'ordre P-1 de la 
prediction lineaire est une fonction d'une stationnarite estimee du signal 
d'entree. 

7. Procede selon Tune quelconque des revendications 1 a 3, dans 
lequel le signal d'entree (x^ est un signal de parole restitue par un decodeur 

(40) a partir d'un flux binaire d'entree (<t>), et dans lequel les parametres de 
prediction du signal d'entree sont extraits du flux binaire d'entree par le 
decodeur. 

8. Procede selon la revendication 7, dans lequel le decodeur (40) 
fournit le residu de prediction (U t ) pour I'adaptation des L coefficients du filtre 

^identification (18). 

9. Dispositif d'identification adaptative d'un systeme (10) auquel est 
applique un signal d'entree, comprenant : 

- une premiere entree (El) recevant le signal d'entree ; 

- une seconde entree (E2) recevant un signal d'observation dont une 
composante est une reponse du systeme au signal d'entree ; 

- un filtre d' identification (18) a reponse impulsionnelle finie representative 
de la reponse impulsionnelle du systeme ; 

- un soustracteur (20) produisant un signal d'erreur en soustrayant du 
signal d'observation le signal d'entree filtre par le filtre d'identification ; et 

- des moyens (22) d'adaptation des coefficients du filtre d'identification en 
fonction du signal d'entree et du signal d'erreur, 

dans lequel le signal d'erreur e t produit par le soustracteur a un 

instant t a pour expression e t = y t -X^H t-1 , ou y t designe la valeur du signal 
d'observation a I'instant t, H M est un vecteur-colonne compose de L 
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coefficients du filtre ^identification, etX^ = (x tl x t _i l --- I x t _ L+1 ) est un vecteur- 
ligne compose des valeurs x t , x t . 1 ,... J x t . L+1 du signal d'entree a I'instant t et 

aux L-1 instants precedents, caracterise en ce qu'il comprend en outre des 
moyens (32, 34, 36) pour obtenir des parametres de prediction du signal 
5 d'entree, minimisant I'energie d'un residu de, prediction sur des trames 
successives du signal d'entree, et en ce que les moyens d'adaptation (22) sont 
agences pour mettre a jour les L coefficients du filtre d'identification (18) en 
ajoutant au vecteur-colonne H t-1 un vecteur-colonne proportionnel a 
e+ 

— = — U t , ou U t est un vecteur-colonne compose des L valeurs du residu 

x t T u t + x T 

10 de prediction du signal d'entree a I'instant t et aux L-1 instants precedents, et X 
un coefficient positif ou nul. 

10. Dispositif selon la revendication 9, dans lequel lesdites trames du 
signal d'entree ont une duree d'au moins 5 millisecondes. 

11. Dispositif seion la revendication 9 ou 10, dans lequel lesdites trames 
15 du signal d'entree presentent des recouvrements mutuels. 

12. Dispositif selon Tune quelconque des revendications 9 a 11, dans 
lequel les moyens pour obtenir les parametres de prediction du signal d'entree 
comprennent des moyens (32, 34) d'analyse par prediction lineaire d'ordre P-1 
du signal d'entree, et dans lequel le vecteur-colonne U t a pour expression 

P-1 

20 U t = X t + 2]a q X t _ q , ou les a q designent les coefficients de prediction lineaire 
q=i 

produits par les moyens d'analyse. 

13. Dispositif selon la revendication 12, dans lequel I'ordre P-1 de la 
prediction lineaire est au moins egal a 5. 

14. Dispositif selon la revendication 12 ou 13, dans lequel I'ordre P-1 de 
25 la prediction lineaire est une fonction d'une stationnarite estimee du signal 

d'entree. 


15. Dispositif selon Tune quelconque des revendications 9 a 11, dans 
lequel le signal d'entree (Xj) est un signal de parole restitue par un decodeur 

(40) a partir d'un flux binaire d'entree (<P), et dans lequel les parametres de 
prediction du signal d'entree sont extraits du flux binaire d'entree par le 
decodeur. 

16. Dispositif selon la revendication 15, dans lequel le decodeur (40) 
fournit le residu de prediction (U t ) aux moyens (22) d'adaptation des 
coefficients du filtre ^identification (18). 

17. Annuleur d'echo adaptatif pour enlever d'un signal de retour (y t ) des 
composantes d'echo d'un signal direct (Xf), caracterise en ce qu'il comprend un 
dispositif d'identification conforme a Tune quelconque des revendications 9 a 
16 t dont la premiere entree (E1) report le signal direct (x t ) en tant que signal 
d'entree et dont la seconde entree (E2) report le signal de retour (y t ) en tant 
que signal d'observation, le signal d'erreur (e t ) constituant le signal de sortie de 
I'annuleur d'echo. 
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